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We study the S — 1/2 Kondo lattice model which is widely used to describe heavy fermion behavior. In 
conventional treatments of the model the Kondo interaction is decoupled in favour of a hybridization of 
conduction and localized /electrons. However, such an approximation breaks the local gauge symmetry and 
implicates that the local /-occupation is no longer conserved. To avoid these problems, we use in this work an 
alternative approach to the model based on the Projective Renormalization Method (PRM). Thereby, within 
the conduction electron spectral function we identify the lattice Kondo resonance as an almost flat excitation 
near the Fermi surface which is composed of conduction electron creation operators combined with 
localized spin fluctuations. This leads to an alternative description of the Kondo resonance without having 
to resort to an artificial symmetry breaking. 
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n the conventional view of heavy-fermion systems, conduction electrons (c electrons) and localized spins (/ 
electrons) hybridize with each other at low temperatures through the Kondo coupling and form a heavy Fermi 
liquid. Thereby, a large Fermi surface develops in this approach, which is determined by the total number of c 
and/ electrons. This picture of heavy fermions is supported by experimental findings, since a narrow quasiparticle 
band is found in the excitation spectrum close to the Fermi energy 1,2 . From a theoretical point of view, a large 
Fermi surface follows within a factorization approximation of the Kondo coupling. In the Kondo lattice 
Hamiltonian 

(1) 



TLq is the kinetic energy of the conduction electrons and Ti.\ is the Kondo exchange (j K > 0). 

S< = { a <>-P /^}flafw an d s i~ Y2v.is ( ff «j8 / 2) c' ia Cijj represent the/and c electron spins. Combining an/creation 

operator from S ; with a c annihilation operator from s ; , and vice versa, expectation values (fipCij\ and (c]pfi,p^ 

can be formed, which lead to a mean- field model of hybridized/and c electrons. A more rigorous treatment of the 
Kondo lattice is based on the so-called large N expansion. Here, N is the/-spin degeneracy, which is artificially 
driven to infinity 3,4 . In the limit °o the key physics is captured as a mean-field theory, and properties at finite 
N are obtained through an expansion in the small parameter l/N. 

Due to the appearance of non-zero amplitudes (f^pC^f) the mean-field description is related to a broken 
symmetry and a hybridization gap in the electronic spectrum. These artifacts are circumvented within our 
approach. A further outstanding drawback of the large N approach is that the cross-over between the heavy 
Fermi liquid and the local moment physics appears as a sharp phase transition where the 1/Nexpansion becomes 
singular 5 . Moreover, the large N approach can not form two-particle singlets for N > 2, such as Cooper pairs and 
spin-singlets. Therefore, anti-ferromagnetism and superconductivity are consequently absent from the mean- 
field theory 5 . Note that the local /-occupation, n[ = Y2afiafi<r' i s a constant of motion of the Kondo lattice model, 
since rl i commutes with the Kondo Hamiltonian. Therefore, the concept of a large Fermi surface in heavy 
fermions was modified by Oshikawa 6 provided that the system can be described as a Fermi liquid. He showed, 
using rather general arguments that the Luttinger sum rule is fulfilled, when also the completely localized spins 
contribute to the Fermi sea volume as electrons. 

However, it is not obvious that the size of the Fermi surface is large in the Kondo lattice model. In particular, in 
the one-dimensional case there is evidence based on Density Matrix Renormalization Group (DMRG) and 
Quantum Monte Carlo (QMC) calculations that the Fermi surface is small 7 9 . On the other hand, for an infi- 
nite-dimensional hypercubic lattice a Fermi-liquid state with a large Fermi surface was found 10 . In this study 
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Dynamical Mean Field Theory (DMFT) and QMC was used avoid- 
ing breaking of the gauge symmetry. Evidence of a discontinuity 
was found in the momentum distribution indicating Fermi liquid 
behavior. 

The DMFT is a widely-used approach to the Kondo lattice which 
relies on an expansion in powers of 1/d, where d is a variable dimen- 
sion 1113 . The idea of the DMFT is to reduce the lattice problem to the 
physics of a single magnetic ion embedded within a self-consistently 
determined effective medium. The qualitative physics of the Kondo 
lattice, including the development of coherence at low temperatures, 
can well be described by this approach. 

Results 

Perturbation theory. The aim of the present work is to provide an 
alternative description of the transformation from localized moment 
to heavy Fermi liquid behavior as a smooth crossover instead of a 
sharp phase transition. To avoid technical details, at the beginning 
let us explain the main idea of our study on the basis of perturbation 
theory. These considerations will be used afterwards to construct 
an appropriate many-particle approach to numerically evaluate 
the electronic spectrum. The one-particle spectral function 
A(k, co) of conduction electrons is defined by A(k,ft>) = (1/2tc) 

e""'dt, where the expectation value and the 



J-oo([M')4< 



time dependence are governed by Hamiltonian (1). The presence 
of the Kondo exchange prevents a straightforward evaluation of 
A(k, co). For that reason we transform the Hamiltonian into a 
diagonal (or at least quasi-diagonal) form by applying a unitary 
transformation to Ti, Ti = e K Ti e~ x . Here the generator X = — X 1 ' 
in lowest order perturbation theory with respect to j K is given by 

k 



1 

X= -j= > 



'£k-£k' 



(2) 



(N L number of lattice sites). The operator quantity /Ckk' = (1 /2) 
X!«/iSk'-k'<rc < ^c£ a Ck7j = (/Ck'k)^ is taken over from the decompo- 
sition of Tti = (j/\/Ni) kk'^kk'- Evaluating the transformation 
to second order in j K , one finds 

Ti-J2 lkC t" Ck " + J2~ Ji i s "' s ) + ^ ( 3 ) 

ka <ij> 

where g k is somewhat changed compared to £k- The second term in 
(3) is the well-known RKKY interaction and E is an additional energy 
constant. Note that in Ti the conduction electrons and the localized 
spins act as independent subsystems which no longer interact with 
each other. After an approximate diagonalization of the RKKY 
interaction, for instance by introducing Schwinger bosons, all 
expectation values formed with Ti can be determined. However, 
the evaluation of A(k, co) also requires the transformation of the 
operator from which the expectation value is taken. This follows 

from the general property (A) = (TiAe^^^j j \Tre~^J = 

(ji^ for any operator variable A, where A = Ae~ x . Thus, we 
obtain 



A(k,w) = ~ 

Z7I 



Ck<r(0>C] 



dt, 



(4) 



where now the expectation value and the time dependence of c ka (t) 
are governed by Tt. An explicit evaluation of c k = e^c\ 



L kr7 



e x up to 



second order in X leads for A(k, co) to the expression 

A(k,oj) = \u k \ 2 8(oj-£ k )+ t^-JZ l v k'k| 2 (Sk-k'-Sk'-k}H,<K'»-£k')> 

(5) 

where the coefficients u k and v k ' k depend on the initial dispersion of 
conduction electrons e k and the Kondo coupling j K . Following 



Doniach's picture of competing energy scales of anti-ferromag- 
netism and heavy fermion physics in 4/ systems 14,15 , in Eq. (5) we 
have neglected the excitations of the localized spin system 
Tij = q/qSq-S_q. Note that in the heavy-fermion regime the 
Kondo scale T K can be considered as large compared to the 
excitation energies of the /system. 

Thus, the one-particle spectrum (5) is built up by a coherent 
excitation with energy e k and amplitude |w k | 2 and by additional exci- 
tations at frequencies e k - with amplitudes |v k ' k | 2 . Usually, such exci- 
tations are responsible for a typical broadening of the coherent 
excitation at E k and determine the lifetime of the quasiparticle in 
the metallic state. In the present case these are the contributions to 
I v k'k| 2 of order j^.. Beyond this second order perturbation theory we 
find that the next higher order oc f K of the incoherent part of Eq. (5) 
leads to the well-known Kondo resonance which has been found by 
Kondo in case of one magnetic impurity. Namely, taking together in 
the second term the parts in j\ and f K Eq. (5) can be rewritten as 



A(k,co) = 8(co — e k ) + 



1 r(co) 



n(co-B k y 



with 



*S(S+l)f K 
4N L 



E 



i+ 



2Jk 



E 



/(«q 



8k' -8, 



q 



<5(eo-£ k ')- 



(6) 



(7) 



where |u k | 2 was approximated by 1. Note that in third order per- 
turbation theory the expectation value (Sk_k''Sk' - k )-f [j ~S(S+ 1) 
becomes local, which means that all interaction effects between loca- 
lized spins at different sites do not contribute since they are of higher 
order. Furthermore note that the sum over q leads to a divergence for 
£ k - = 0, and the additional excitations in A(k, to) become dominant. 
They are located close to the chemical potential co = 0. Since their 
energy is momentum- independent, i.e. independent of k, this set of 
excitations has to be interpreted as the well-known Kondo res- 
onance. The approximate result (7) agrees with the imaginary part 
of the Green's function for the impurity Kondo model in perturba- 
tion theory (compare e.g. with Ref. 16), where A(k, w) = (l/n) 
9G(k, w). Using moreover the Green's function G(k, co) in the gen- 
eral form 

G(k,a»)= — (s) 

co — Ek" iT(co) 

(real part of the self energy is neglected) one also finds that F(co) has 
to be interpreted as imaginary part of the one-particle self-energy. 
Thereby, higher order terms in the denominator of A(k, w) = (l/ri) 
9G(k, w) are neglected 17 . 

Note that the perturbation expansion (7) of T(m) in j K breaks 
down below a characteristic temperature when the second term in 
the bracket becomes of the order of the first term. This temperature is 
often used to define the Kondo temperature and is conventionally 
obtained from an expansion of an effective Kondo coupling 18 . To see 
this the sum over q can easily be evaluated 



r(co) 



r(co)- 



nS(S+l)f K 



nS(S+l)f K p 0 



max(|co|,T) 



D 



-o(f K ), (io) 



where in the last line also the integration over k' was performed. Here 
p 0 is the electronic density of states. In a final step the single-impurity 
Kondo temperature T K = D exp( — 1/jgPo) i s introduced leading to 



I» = 



nS(S+l)f K p 0 



3-2pa/Km 



max(|ft)|,T) 
Tk 



(if) 
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One should again emphasize that no difference from the single 
impurity case 16,17 has occurred up to now. 

The previous considerations are based on perturbation theory 
leading to a divergence in A(k, co). This divergence can be removed 
by including higher order contributions to infinite order inj K . For the 
single impurity model, higher order perturbation theory leads to 
equally divergent contributions. This series was summed up by 
Abrikosov 19 in the leading ln( T/T K ) approximation to give the fol- 
lowing result for the T-matrix, T(e)~(j K /Ni)(l — (j K p 0 /N 
In(|«|/D)))-Vs. 

Non-perturbative many-particle approach to the 2d-Kondo 
lattice. The previous perturbation treatment leads to two 
characteristic problems in the evaluation of the one-particle 
spectral function A(k, co). (i) There are contributions to the 
spectral intensity which diverge and lead to unphysical results, (ii) 
There is no difference from the single impurity case in lowest order in 
)k- These problems can only be solved by including higher order 
contributions to infinite order in j K . In the following, we therefore 
apply the Projector-based Renormalization Method (PRM) 20,21 to 
evaluate A(k, co) for the Kondo lattice. The general concept of the 
PRM (see methods section) is similar to what has been discussed 
above: To enable the solution of Hamiltonian (1), the Kondo 
interaction is again integrated out. However, a sequence of unitary 
transformations is used instead of a single transformation as in 
perturbation theory. Thereby, the PRM ensures a well-controlled 
disentanglement of higher order interaction terms leading to an 
effective Hamiltonian of the same operator structure as in Eq. (2). 
However, the renormalized parameters £ k , /« and E are determined 
self-consistently taking into account contributions up to infinite 
order in the Kondo coupling. Moreover, within our non- 
perturbative approach expression (4) for A(k, co) is still valid. The 
amplitudes u k and v k * k are also renormalized under the influence of 
higher order contributions and will be marked by tilde symbols as 
well. At this point we would like to emphasize the main motivation to 
develop our own theoretical treatment: Firstly, we note that the 
perturbative result (6) is also accessible within the PRM by 
expanding the renormalization equations up to lowest order. 
Furthermore, our approach has important advances over other 
theories for the Kondo lattice: The PRM neither suffers from 



additional constraints, as in the large N approach, nor from an 
expansion for large dimension d as in the DMFT. Thereby, the 
local /-electron occupation remains a constant of motion in the 
calculation. 

In the actual evaluation of the renormalization equations we 
restrict ourselves to two dimensions in order to minimize the numer- 
ical effort. Also, we consider electron concentrations only away from 
half-filling which is a special case since there the system becomes 
insulating. We have solved the renormalization equations self-con- 
sistently on a square lattice with N L = 10 3 lattice points. Fig. 1(a) 
shows A(k, co) for T = 0 in a very small k-region around the Fermi 
momentum k F for electron filling n c = 0.3. Clearly one can see the 
coherent excitation e k which crosses the Fermi level at co = 0. Note 
that the renormalization of e k is very weak. This behavior is found to 
be always true away from half-filling. In contrast, at half-filling n c = 
1/2 renormalization contributions to e k become strong due to nesting 
effects of the Fermi surface and lead to the opening of a gap at the 
Fermi energy. This results in an insulating phase 22,23 and will be 
discussed in a forthcoming publication. Moreover, as expected, an 
additional almost k-independent excitation is found at the Fermi 
level which corresponds to additional contributions described by 
the second term in Eq. (5). In contrast to the perturbative treatment 
the renormalized amplitude vvk is no longer divergent in the PRM. 
However, it is still dominant for k' ~ k F and any value of k. As before, 
the multitude of these excitations has to be interpreted as Kondo 
resonance. It results from the second term in A(k, co) and not from 
a hybridization effect off and c electrons. This important feature of 
the Kondo resonance can also be recognized in Fig. 1(a). In contrast 
to the coherent excitation, for the Kondo resonance we find a con- 
tinuous spread of spectral intensity at the Fermi level throughout the 
whole momentum region. Due to the finite number of lattice points, 
considered in our calculation, the one-particle kinetic energy e k is 
defined only for a finite number of k-points. Therefore the coherent 
part of A(k, co), following the dispersion according to the first term of 
Eq. (5), shows a set of discrete bright spots (19 visible in Fig. 1) which 
is characteristic for a finite system. In contrast, the second part in Eq. 
(5) which contributes to the Kondo resonance is made up of an 
internal sum over all momenta k' leading to the observed continuous 
momentum character. Considering finite electron systems, this char- 
acteristic feature could possibly be picked up by experiments to find 



(a)T = 0 (b)k B T/(10" 3 t)=5 




-0.8 -0.4 0 0.4 0.8 -0.8 -0.4 0 0.4 0.8 

(k - k F ) / (1 0" 2 ti) (k-k F )/(10 _2 7t) 

Figure 1 | (a) One-particle spectral function A(k, co) for n c = 0.3 and j K /t= 0.3 as afunction of co for fixed values of k = k(e x + e y ) in a small wave-vector 
region around k F . The k-resolved Kondo resonance peak around co = 0 is clearly seen for all k values, (b) A(k, co) for the raised temperature k B T = 5 • 
10~ 3 f, which is above the Kondo temperature T K . For fixed k around k F only the coherent excitations survive whereas the Kondo resonances disappear. 
Note that the discrete bright spots for the coherent excitations are characteristic for finite systems which is not the case for the Kondo resonance. 
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experimental evidence for the different nature of the Kondo res- 
onance. In Fig. 1(b) the spectral function is given in the same 
wave-vector region as before for the raised temperature of T = 5 • 
10~ 3 t. Now the Kondo resonance has disappeared and only the 
coherent excitation has survived, which means that T is higher than 
the Kondo temperature T K . Let us mention that the overall behavior 
in Fig. 1 looks similar to numerical QMC results 24,25 for the periodic 
Anderson model in one dimension for large values of the Coulomb 
repulsion U. As is well known, the Anderson model reduces to the 
Kondo model in the large U limit and low lying energy level ty of the 
localized/electrons. In particular, in Ref. 24 some /-like bands with 
low weight were found very close to the chemical potential, two of 
which seemed almost not to disperse which is quite similar to the 
behavior of the resonance mode in Fig. 1. 

Fig. 2 shows the density of states p c (co) = (l/N L )2 k A(k, co) as a 
function of co for three different values of T in a small frequency 
region around co = 0. For the two lower Tthe Kondo peak around co 
= 0 appears as a consequence of the Kondo contributions in A(k, co). 
The extra background in the spectrum is mainly composed of coher- 
ent excitations with respect to different momenta k. Additional inco- 
herent excitations with energies different from zero may contribute 
to the background as well. We have also checked the sum rule 

A(k,co)dco = 1, which at first glance seems to be contradicted 

00 

by the obvious T-dependence of the different values in Fig. 2. 
However, the frequency interval outside the restricted cu-region of 
Fig. 2 also contribute to the sum rule so that the sum rule is almost 
perfectly fulfilled for the different T curves. Similarly, particle-hole 
symmetry seems to be present for the curves of Fig. 2. However this is 
not the case for n c = 0.3 when the full frequency region is considered 
(not shown). 



Discussion 

Firstly, let us discuss how the Kondo scale enters the formalism. To 
address this point, we start from relation 5Zi(Sf s i) = 
5Zkk'(^ k ' k )' wn i cn suggests that the expectation value (/Ck'k) parti- 
cipates in the formation of the singlet state at low temperatures. In 
Fig. 3 the expectation value (/Ck'k) is shown for wave vectors close to 
the Fermi surface, k' = k ~ k F . As expected, (/Ck'k) turns out to be 
negative. Moreover, it exhibits a strong increase in magnitude with 
decreasing temperatures. This behavior is an indication for the 



100 
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60 
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— k B T/(10 t) = 5 

— k B T/(10' 3 t) = 0.2 

— T = 0 




(0/(10 f) 



Figure 2 | Density of states p c (a>) for n c = 0.3 and j K lt = 0.3 as a function 
of oj for three different temperatures, plotted for a small frequency 
region around co = 0. Clearly seen is the Kondo resonance at oj = 0 for the 
two lower temperature values T = 0 and T = 2 ■ 10~ 4 t. For T = 5 • 10~ 3 t 
the resonance peak disappears since T is larger than the Kondo 
temperature T K ~ 5 • 10~" t (see text). 



formation of the singlet state, to which predominantly wave vectors 
k', k close to the Fermi surface participate. In real space this gives rise 
to correlations (S ; • Si) between local and conduction electron spins 
which are limited to a range of about 50 lattice sites (see inset). For 
comparison, Fig. 3 also shows the complete expectation value (S, ■ s,-) 
as a function of T. Obviously, (S, ■ s,) is negative as well but almost 
independent of T. Because of the former relation between (S, ■ s,) and 
(/Ck'k) one concludes that expectation values (/Ck'k) f° r wave vectors 
away from k P are not affected for low temperatures. These wave 
vectors form the majority in the Brillouin zone, which explains the 
marginal T-dependence of (S,- • s,). The width of (/Ckk ) in Fig. 3 can 
be identified with the Kondo scale k B T K and is of order 5 • 10~ 4 t. The 
same energy scale is also found in the T-dependence of the electronic 
part of specific heat Cy(T), which can be evaluated in the PRM as well 
(not shown). For low T a strong increase ~ T is found which is 
usually associated to a large heavy fermion mass. For T > T K the 
linear T-slope changes to a much lower value as known from normal 
metals. 

In conclusion, we have applied the Projector-based Renormali- 
zation Method (PRM) to the Kondo lattice model in d = 2 in order 
to study the electronic one-particle spectrum in the heavy-fermion 
regime. As main result of our study, we have found that the Kondo 
resonance is caused by additional low-energetic contributions to the 
spectral function A(k, co), which are dispersionless and distributed 
continuously in momentum space even for a finite electron system. 
Such a continuous distribution disagrees with a mean-field descrip- 
tion of the Kondo resonance and proves that the Kondo resonance is 
a true many-body effect. 

Looking at the electron density n c , one finds that n c is related to the 
electronic spectral function via n c = (I/Nl) J A~ (k,co)dco. On 
the right hand side, all excitations from the non-perturbative exten- 
sion of Eq. (5) contribute. This relation fixes the position of the Fermi 
energy, when n c is fixed. Using our result for A(k, co), the Kondo-type 
excitation turns out to be located slightly below the Fermi energy, i.e. 
inside the Fermi volume, which is consistent with photoemis- 
sion experiments. This low-energy quasi-particle is the famous 
Abrikosov-Suhl resonance indicating heavy fermion behavior which 
is mainly characterized by a large effective mass of conduction elec- 
trons. The resonance emerges from the second term in Eq. (5) in 




1 2 3 4 5 

k B T/(10" 3 t) 

Figure 3 | Spin coupling (/C^k) for k' = k ~ k F as a function of T (black 
curve). For Tbelow T K (/C k 'k) drops to rather small values which indicates 
the formation of the heavy-fermion state. In contrast, (S, ■ s,) is almost T 
independent (red curve). This behavior demonstrates that only wave 
vectors close to k F contribute to the formation of the heavy-fermion state. 
The inset shows the correlation function (S, • Sj) as function of (i — j). 
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addition to the conventional coherent contribution. Thus, experi- 
mentalists might interpret both excitation types, i.e. the coherent 
excitation £k and the Kondo resonance, as ingredients of the large 
Fermi volume within a Fermi liquid picture resulting from a mean- 
field description. 

Unfortunately, the question whether the Kondo lattice in 2d is a 
Fermi liquid can not be answered by our present treatment. On the 
one hand, for the electron filling n c = 0.3 considered in this work the 
results are comparable with the Kondo impurity model, which is 
known to be a Fermi liquid 19,26,27 . The Kondo resonance always arises 
at £k ~ 0, i.e. at w ~ 0 and is wave-vector independent as in the single 
impurity case. Moreover, we have found that the wave-vector 
dependence of the spin correlation function (Sk_k' ' Sk-_k) is rather 
smooth, i.e. the correlations between local spins at different sites are 
rather weak. Therefore, we conclude that for parameter values not 
too close to the quantum phase transition to antiferromagnetism, the 
impurity physics largely governs the Kondo physics of the lattice 
model at n c = 0.3. On the other hand, one should emphasize that 
correlation effects between different lattice sites were found in the 
spin correlation function (S, ■ S;) which is spatially extended in a 
range of 50 lattice sites (inset of Fig. 3). Thus, long ranging coherence 
effects between local and conduction electron spins do not disturb 
the impurity-like behavior of the Kondo lattice. Moreover, for k 
values on the Fermi surface, k = k F , the numerical evaluation of 
the quasiparticle weight |«k| 2 leads to a very small value. A non- 
vanishing of the quasiparticle weight |«k| 2 (usually called Z k ) would 
be an important evidence that the Kondo lattice in d = 2 has the 
character of a Fermi liquid. 

A possibly simple way to check experimentally the Fermi liquid 
behavior is to measure the k-dependent occupation number. Note 
that at finite temperatures the expected jump at k F for T = 0 is 
smeared out. Therefore, one has to extract the finite temperature 
effects by use of an adequate Fermi distribution function, which is 
the usual way to interpret photoemission experiments at low tem- 
peratures. A jump at k F would support the Fermi liquid picture. 

Methods 

Projector-based renormalization method (PRM). The PRM starts from a 
decomposition of a many-particle Hamiltonian Tt into an unperturbed part Tio and 
into a perturbation Hi. The latter part accounts for transitions between the 
eigenstates of 'Ho ■ The basic idea of the PRM is to integrate out the perturbation H i by 
a series of unitary transformations. For practical applications the unitary 
transformations are best done in small energy steps AX. Thereby, the evaluation in 
each step can be restricted to low orders in "Hi. This procedure usually limits the 
validity of the approach to parameter values of H i which are of the same magnitude as 
those of Hq. Assuming that all transitions with energies larger than some energy 
cutoff X are already integrated out, the transformed Hamiltonian Tix for the Kondo 
lattice model reads, Jix — TCo^ + Hi.x with 

W <U = ^2 ^4*%* + XXA'S-q + B;. (12) 

kff q 

Hy= —}= jkk,;.® (^ - - 8m )K-n, (13) 

V N L kk 

where the operator quantity /Ckk was introduced before. Due to the elimination of 
high-energy transitions above X the coefficients in Eqs. (12), (13) depend on X. As in 
the perturbative treatment, an effective exchange interaction between the/ spins as 
well as an energy constant £; and a new k-dependence of the Kondo coupling;^ are 
generated by the transformation. The ©-function in T~i\j limits the excitations to 
energies smaller than X. To find the /.-dependence of the coefficients in H ; we study a 
subsequent unitary transformation of 7ix to a new Hamiltonian Ti^-Ax 

H ? _- Ai = e x ^H l e- x >~»'-. (14) 

Thereby, all transitions between X and a somewhat reduced cutoff X — AX will be 
eliminated. The generator Xx_aa of the unitary transformation 

X A ,AA = —7^= ^k'k('^AA)/Ck'k 

(15) 

X kl[ ( X, AX) = 7k ' M Qkk,A 1 1 ~ ®k-M -ax), 
ek,a — %a 

(©k'M = ©(^ — \ E k',A ~ £ m|)> i s fixed by the requirement that Jii-AA n ° longer 
contains excitations with energies larger than a — A?.. Since the elimination of 



transitions is confined to a small energy shell AX, the evaluation of TCa-AA from ?i A 
can be restricted to low order contributions in jk'M- Expression (15) is a slight 
generalization of the generator X, which was used to derive 7i in Eq. (2). By evaluating 
transformation (14) to order jj^, x together with (12), (13) we obtain discrete 
renormalization equations (or flow equations) connecting the coefficients at cutoff X 
with those at cutoff ). — AX. For instance, the renormalization equation for jk'k ^ reads 

MU-AA -jk%A = r kq;k'q( 2 ( <: ^ C 1 a )) 1 ( 16 ) 

1 q 

where r;j, q . k , k - ©qq^qq'A'k^A^) + ©kkvjkk',^q'q(MA). The complete 
elimination procedure starts from the original model (1), where the coefficients are 
fixed to e^a — £ k> 7q,A — 0> £a ~ 0>7kk',A — j K . Here X = A denotes the largest 
transition energy of the original Hamiltonian (1). Proceeding insteps A?, until X = Ois 
reached all transition operators from Tti will be used up. One arrives at the final 
Hamiltonian Tt A =a> which has exactly the same form as the Hamiltonian 7i of Eq. (3). 
Thus, from now on quantities with tilde symbols will always denote the fully 
renormalized quantities at cutoff a — 0. Note that the explicit evaluation of the unitary 
transformation leads at first to operator expressions which differ from those of the 
generic ansatz (12), (13). Therefore, an additional factorization approximation has to 
be performed in order to trace back these operators to the generic ones. As a result, the 
renormalization equations also depend on expectation values {■•■}, which have to be 
solved self- consistently together with the renormalization equations of Ha- For the 

evaluation of these expectation values we use the invariance property (_4) — ^_4^ _. 

To find the desired spectral function A(k, to) we start from Eq. (4) with 

c k(7 — c\^ G {). — 0). The A-dependence of c ko .(/i) is calculated via the following ansatz: 

C L W = U^Atr + -T^T Sk k ' ~T ^ ( 17 ) 

with initial conditions w^a — 1, vvk,A ~ 0- The operator structure of (17) is suggested 
from the first order expansion in Xx,aa- The renormalization equations for the X- 
dependent coefficients u^x ana " v k'k,;. are obtained by transforming c\. a (?„) to 
c k(T (A — AX) in analogy to the transformation of TC^- The equation for Vkjc',/1 to first 
order in X^aa reads 

Vk',M-A,l - Vk',M ^ U KX X k ' k (lA?.) + V qk,A'q(MA) ^C^C^ l) , (18) 

The renormalization equation for can be evaluated in the same way. Alternatively, 
it can be found from the sum rule / j^a-(A),c£ a (X)^ \ — 1, with I replaced by a — 
AX 

\Uk,A-AA\ 2 ^ 1 - ^fYl | V qM-AA| 2 {Sk-q-S q _k) ~ 
L q 

l (19) 

— ^XS k ^- AAV q'k.A-AA<£qq')- 
*- N L qq' 

Note that due to (18) the coefficients v^x-aa ana " v q\k,A-AA in (19) at cutoff X — AX 
can be expressed by Vk',M and u^ - Thus, equations (18), (19) relate the two quantities 
"M-aa an d v k'XA-AA at cutoff X — AX with the same quantities at X, Since the 
renormalized Hamiltonian Tix=a — 7i has the same decoupled form (3) as in 
perturbation theory, the renormalization equations (18), (19) together with (17) lead 
to the final result for the spectral function A(k, m). Thereby, A(k, to) takes the same 
form as (5) but with the fully renormalized parameters of Hamiltonian H and of c kff . 

As mentioned above, the Kondo resonance results from additional contributions to 
the spectral function A(k, tu) close to the Fermi surface. Looking at expression (5) for 
A(k, to) one notices that only terms from the sum over k' can contribute to the Kondo 
resonance which have k' ~ k/?, so that £k ~£k F =0. At the same time, also the pre- 
factors |vk',k| 2 become large for the same k' values. This latter property can be 
understood by means of renormalization equation (17) for v k 'k,A : Taking k' ~ k^and 
arbitrary values of k kj.-, the first term u^X^kiX, AX) on the right hand side of ( 1 8) is 
nonzero but small since X^-^X, AX) is small. This follows from the structure of Xfk£/l> 
AX) in (15). Similarly, in the second term of Eq. (18) those terms from the sum over q 
contribute most for which q ~ kp, since then Xk' q ( X, AX) is large. Thus, one concludes 
that terms with wave vectors k' ~ k F lead to a finite value of Vk'ku in the renorma- 
lization procedure for X — > 0. This property explains the appearance of the Kondo 
resonance in A(k, to), which is almost independent of k [Fig. 1]. On the other hand, 
from sum rule (19) one can show that for X — > 0 the renormalized amplitude |iik| 2 of 
the coherent excitation in A(k, m) becomes reduced. For the case k — kp, which was 
excluded above, the numerical evaluation of the renormalization equations (18), (19) 
leads to a very small value of the coherent coefficient |wk|"> which can not be dis- 
tinguished from zero within the numerical accuracy. 
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